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C/3 

^^ ' It has hitherto been widely considered that a mixing layer is unstable at all Reynolds 

j^^' numbers. However this is untenable from energy considerations, which demand that there 

must exist a iron-zero Reynolds number below which disturbances cannot extract energy 

s^ from the mean flow. It is shown here that a linear stability analysis of similarity solutions 

f^ ' of the plane mixing layer, including the effects of flow n on-parallelism, using the minimal 

o' 

\l ' composite theory and the properties of adjoints following 



Govindaraian fc Narasimha (2005|) . 



o 
\o 

^— N ' resolves the issue by yielding non-zero critical Reynolds numbers for coflowing streams of 

(\ ' any velocity ratio. The critical Reynolds number so found, based on the vorticity thick- 

• 1— I ' 

J>^' ness and the velocity differential as scales, varies in the narrow range of 59 to 55 as the 

^' 
Mh' velocity ratio goes from zero to unity. 

> 

X 

1. Introduction 

The motiv ation for the present work ar ises from an analysis of the stability of a mixing 



layer due to 



Betchov fc Szewczvk ('1963I) . This analysis, based on the Orr-Sommerfeld 
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(OS) equation, showed that the neutral curve in the wave-number (a) - Reynolds number 

{Re) plane approaches the origin a = a,s Re goes to 0, leading the authors to conclude 

"No minimum [i.e. critical] Reynolds number is fou nd" for the flow. Mixin g layers have 



been studied extensive ly since, a nd recent texts, e.g. 



Criminale et al. (20031 p. 90) and 



Huerre fc Rossi (19981 pp. 81-294) 



Drazin fc Reid (2004 p. 199-201), stiU report the crit- 



ical Reynolds numb er as zero. Th is is intrig uing, for several early studies of stability 



(e-g- 



Prandt] (198,^ pp 178-183), 



Lin fl95,'Tl pp. 31-32, 59-I-), going back to the two- 



dimensional analysis of lLorentz (19071 pp. 43-71)), show that (under certain reasonable 
conditions) a two-dimensional incompressible viscous flow must be stable at sufficiently 
low Reynolds numbers. Astonishingly, however, there is no analysis in the extensive lit- 
erature on hydrodynamic stability that yields a non-zero value for this critical Reynolds 
number. The mixing layer being a very basic flow type, the absence of any definitive 
commentar y on the above issue has led other studies of flows modelled on the mixing 



layer (e.g. 



Solomon. Hollowav fc Swinnev 1993: 



Ostekin. Cumbo fc Liakopoulos 1999|) 



to base their analyses on the assumption that th e critical Re yn olds number is zero. 



Esch (19571) to 



Balsa (1987|) . start with 



Now all analyses of mixing layer stability, from 
the OS equation, which is valid only for strictly parallel flow. However, the width of a 
laminar mixing layer scales as x^'^ (where x is the streamwise distance), so the rate of 
change of thickness (and hence also any parameter measuring the degree of flow non- 
parallelism) becomes infinite in the limit of x (equivalently Re) going to 0. In other 
words, existing studies are based on the assumption of no non-parallelism in a situation 
where any measure of non-parallelism would be infinite. Thus flow non-parallelism may 
be expected t o play a crucial role in determining the stability characteristics in the 



limit. Though 



Bun fc Criminale (1994(1 note that "with viscous effects, the basic flow 



should be treated as non-parallel", and other texts (e.g. 



Drazin fc Reid 20041 p.l97) 



Stability of a laminar mixing layer 3 

emphasize that "at smah values of Reynolds number the parallel-flow assumption is of 

questionable validity" , no investigation accounts for the non-parallelism. We show here 

that a consistent non-parallel flow theory yields finite, non-zero critical Reynolds numbers 

based on vorticity thickness and velocity differential in the range 55 to 59 depending on 

the velocity ratio. 

In a departure from earlier work, the similarity solution of the laminar mixing layer is 

used as the base flow everywhere in the present analysis. The earlier studies, beginning 

with ones that assumed zero viscosity, have moved from th e discontinuous profile due t o 



Helmholtz to the hyperbolic tangent profile c onsidered by 
Examples of other approximations appear in 



Betchov &: Szewczvk Cl963|) . 



Esch ('I957f) . and all later studies to our 



knowledge use one of these. The similarity profile is a better approximation of reality 
than any of the above, and its use through what we have called minimal composite theory 
is both appropriate and convenient in the present approach. 

The rest of the paper is arranged as follows. In §|21the similarity solution of the mean 
flow profile is presented and some remarks follow. The stability problem is posed and 
the method of soluti on briefly outlined in jjEl essentially follo wing the approach used for 



boundary layers by 



Govindaraian & Narasimha fl997l 20051) . In §01 the results of the 



stability analysis are presented and compared with earlier work. 

2. Similarity solution of the basic flow 

A plane incompressible mixing layer that develops in the positive x-direction is consid- 
ered. The two free-streams flow with velocities Uoo and XUoo (with A < I) before coming 
into contact with each other at the origin {x = 0). Both streams are semi-infinite in their 
lateral (y) extent. We omit the case of counterflowing streams (A < 0) from the present 
discussion, since they are not relevant to the non-parallel flow theories. Since no other 
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external length scale is present in the problem as formulated, after a sufficient distance 



from the origin similarity may be taken to be valid. 

Henceforth all variables subscripted by d are dimensional quantities. The streamwise 
coordinate and the coordinate in the direction normal to the flow are Xd and y^ re- 
spectively, and v is the kin ematic viscosity. The similarity analysis here differs from 



Schlichting fc Gersten f2fl04L pp. 175-176) only in the definition of the length scale Id in 



the direction normal to the flow, 



ld{xd) = \jjj^- 



(2.1) 



The streamfunction is defined as 



^d{xd,Vd) = Uookixd)^{y), 

y — Vd/ld being the similarity coordinate. The momentum equation in terms of the 
non-dimensional streamfunction <P becomes 



$'" + }.,p<p" = 0, 



(2.2) 



where the primes denote derivatives taken with respect to y. The boundary conditions 
to be satisfied are 

$'{y = +oo) = 1, 



<^'{y = -oo) = A and 



iy'^'-i')]y=oo + ^iyi''~ny=^oo 



0. 



(2.3) 



The third boundary condition has been a subject of some controversy. As formulated 
above, it is derived from matching the pressure across the mixi ng region, following 



Ting fl959|) . The zero net transverse force con dition suggested by 



von Karman fl92l|) 



results in an identical third boundary condition. iKlemp fc Acrivos (1972(1 pointed out an 
inconsistency in the above formulation and showed that this condition is still incomplete, 
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Figure 1. Stream-wise velocity profiles from the solution of 12.211 for diflerent velocity ratios. 
Solid line, similarity profile for A = 0; circles, i(l+tanh{2(y+l)/5}); dash-dotted fine, similarity 
profile for A = 0.5; short dashes, similarity profile for A = 0.75. 

and further that it cannot be resolved within the context of classical boundary layer 
theory. Using a different ap proach to formulate asy mptotic expansions to solutions of 



the Navier-Stokes equations, 



Alston fc Cohen fl992ri have however argued that this con- 



dition remains the best option compared to the other alternatives used in the literature. 

The velocity profiles obtained for different values of A are shown in figure^] Apart from 

the displacement of the dividing strea mline in the negative y direction, these solutions 

are identical to those obtained first bv iLoc 



A = 0, was validated against that given in 



t (19511). The result for the 



half jet, i.e. for 



Schlichting fc Gersten ('2004 pp. 175-176) 



A suitably transformed hyperbolic tangent function is superimposed on the similarity 
profile for A = for comparison. 

Note that for every flow configuration with a A in the range to 1, there corresponds 
an identical but vertically flipped flow with velocity ratio 1/A. So it is expected that the 
similarity solution becomes increasingly symmetric about the x axis as A approaches 1. 
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In the limit A — > 1, redefining the non-dinicnsional streamfunction as 

<f(y) = (l-A)/(y)+y, 

the similarity equation H2.2|l becomes 

/"' + ^[y + (l-A)/]/" = 0, (2.4) 

with appropriately transformed boundary conditions. Writing (1— A) as a small parameter 
e and assuming a solution in the form of an asymptotic scries in powers of e, 

/ = /o + e/i + e^/2 + (higher order terms), 

we get the streamwise velocity correct to second order in e as 

Udoixd,yd) = Uoo[l + e./^(y)] = C/oo|l - K^l - CTiiy/2)) \ + 0{e^). 

While such analytically expressible profiles arc ubiquitous in the literature, our purpose 
is to demonstrate that the similarity solution smoothly merges into the error-function 
profile for A close to 1. It is therefore not expected (as we shall confirm shortly) that 
the stability characteristics change very much between these choices. This puts into 
perspective one of the nuances of our analysis that is different from existing ones, for it 
must be emphasized that the marked deviation observed in the final result is not to be 
attributed to the differences in the assumed velocity profiles but to the non-parallel flow 
analysis. 

3. The non-parallel stability problem 

A brief outline of the minim al composite theorv is described in w hat follows. A review 



of the method can be found in JNarasimha fc Govindaraian (^200C ) 



Each flow quantity, for example the streamfunction, is expressed as the sum of a mean 
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<Pd and a perturbation (pd^ where 

^d{xd,yd) = Uoold{xd)<^{y) and (3.1) 



(i>d{xd,yd) = Uoold{xd)<P{x,y)exp 



i / a{x)dx — Lut 



(3.2) 



with the (complex) phase speed of the disturbance c = Lu/a. Inserting these into the 
Navier-Stokes equations for two-dimensional incompressible flow written in terms of 
the streamfunction, and retaining all terms nominally upto 0{Re^ ), the non-parallel 
stability equation can be written as 

Af{(j)} = 0, with boundary conditions 0, Di^— »0 ai y — ±00, (3-3) 

where, for the mixing layer, 

AA = (^' - c)(D2 - a^) - $'" + 4-/d4 - 2a^D^ + a^ + p^T)"" 

ane y 

+p(<Z>'D2+ <Z>"D) -a^[2py{^' -c)+p$]B-pa'^c + p^"' 

+ {i<P' - cjReaa' + \<P"' + a^(3$' - 2c) - -Z-'D^li^e— I, 

ox } 

D is the derivative with respect to y and p = Re{dld/dxd) — 1/2. The Reynolds number 

is based on Uoo and Id- Note that in the minimal composite theory the operator TV is flow- 

specific, and so is different e.g. fro m that for the boundary layer (see below). Following 



Govindaraian fc Narasimha fl997J) the non-parallel operator TV is expressed as the sum 
of an operator M that contains all the lowest order terms and an operator 7i comprising 
the higher order terms, 

ane 

n = — \p<PD^ - 2a^D^ + a"^ + pi'P'D'^ + <P"D) - a^\2py{<P' ~ c) + p<P]D 
Re 



-pa^c + p'^"'\ + {i<P' ~ c)aa + S 
where S = [<?'" -I- a'^{3<P' - 2c) - <P'V>% 



d_ 
dx' 
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The relative order of magnitude of each individual term varies with y because of the 

presence of the critical layer at <?' = c. While constructing the minimal composite 

equation that yields results correct upto 0{Re^ ), the order of magnitude of any term 

within the y-domain is considered. The final equation consists of all terms that are 

at least 0{Re^ ) somewhere, and rejects all that are o{Re~ ) everywhere in the do- 

main. It must be emphasized th at the above equations are different from those in 



Govindaraian fc Narasimha f 2001111 . which were formulated for a boundary layer. For 
instance, in the construction of the operator A4 for the present case clearly no wall layer 
considerations are needed. Furthermore, the termp<?D^, being of higher order everywhere 
in the mixing layer, is in the operator 7i, whereas in the boundary layer it is part of A4. 
The solution procedure for estimating the growth of the disturbance, though, remains 
the same, and is given in that paper in detail. For reference, we recapitulate some of its 
essential points below. 

It is instructive to note that though the non-parallel operator has partial derivatives 
in both X and y, the lowest order terms (comprising M) contain derivatives only in y. 
Taking this as a cue, the total solution (f> is expressed in two parts, 

(^(x, y) = A{x)(t)„i{x, y) + e(l)h{x, y). (3.4) 

Here A(x) is the amplitude function that captures the streamwise variation of the lowest 
order solution 0m, which satisfies the equation 

M{(j}„^] = 0. (3.5) 

We consider the downstream growth of disturbances at a constant value of the similarity 
variable y. For future reference, we also define a complex effective wavenumber associated 
with the non-dimensional streamwise disturbance velocity u as 

yid^ , 1 d{d<j,^/dy) \ 
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such that the growth rate g of -u is given by 

g = -Im(a.ff) ^ -a. + Rc(^-- + g^^Jg^^^^ ) - (3-7) 

where ai is the imaginary part of the full complex wavenuniber a. We find that, in general, 
a disturbance may amplify at one y and decay at another. Moreover, one disturbance 
quantity could be amplifying while others decay. The derivative of 0m with respect to x is 
obtained by solving (|3.5|l for a nearby Re and noting that d/dx — pd/dRe. Substituting 
from (|3.4|) into (|3.3|) . and noting that Af — M +H, we obtain the amplitude evolution 

equation 

dA 

AH<i),^ + —S(j)r,^^-eM(j)h, (3.8) 

the truncated terms here being 0{Re^^'^) compared to the largest of the retained terms. 
Further the expression for the adjoint of the operator M is found to be 

M = (^' - c*)(d2 - a*2) + 2<Z'"D ^D"* + 0(Re~^/^), 

^ '^ ' a*Re 

where a quantity superscripted with an asterisk represents its com plex conjugate. Us- 



ing the property of adjoints (cf. 



Govindaraian fc Narasimha 2005|) . the contribution to 



the growth g due to change in the amplitude function, i.e. the quantity (l/yl)dyl/da;, is 
calculated without the need to specifically compute the higher order solution 0/i. Equa- 
tion (|3.8|) therefore reduces to 

/oo dj4 (°^ 

X*'H{<l^rn}dy + — x*'5{0m}dy= (higher order terms), 

-oo ^-^ J ~oo 

where x is the solution of the adjoint problem. 

Mix] = 0. 

Now A ^ 0(1) and dA/dx is 0(iZe^^), so the error in the estimation of the growth 
rate will be o{Re^ ). But as we integrate over a large streamwise distance of 0{Re), the 
error in the amplitude is expected to be o(l). Here, a solution of the partial differential 
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equation H3.3|l to the desired order of accuracy has been obtained by solving a parametric 

ordinary differential equation and using the property of adjoints. An obvious advantage 

is that the complexity of the problem is significantly reduced when compared to the 

exercise of solving it as a full partial differential equation. 

4. Results 

4.1. Accuracy of results 

Both the parallel and non-parallel stability analyses involve solving an eigenvalue prob- 
lem. The conditions at infinity are posed at distances sufficiently far away to ensure that 
results are independent of the size of the domain. For example, on the curve of marginal 
stability, the Reynolds number for a given a is required to be identical up to the fifth 
decimal place before any further increase in domain size is considered unnecessary. We 
notice that at large y, <P"' vanishes and the higher derivatives of the eigenfunction with re- 
spect to y (appearing in the viscous term) become smaller. Therefore, as an estimate, the 
eigenfunction decays as exp(±ay), i.e. its e-folding rate depends on the wavenumber of 
the disturbance. The y-domain thus needs to be larger for smaller wavenumbers. Further, 
since the eigenfunctions are discretized as eigenvectors, a sufficient number of grid-points 
must be contained in the domain so that the results are independent of resolution. 

Now, an increase in the number of grid-points directly increases the size of the matrices 
involved in the eigenvalue problem and this affects the computational effort adversely. For 
non-dimensional wavenumbers larger than 0.2, we impose outer boundary conditions at 
y — ±20 and use 81 grid-points with a sinh grid-stretching. For the smallest wavenumbers 
considered {a ~ 0.02) the conditions at infinity are imposed 'At y — ±300 with 150 grid- 
points. Intermediate values were taken to expedite the solution process whenever this 
did not affect the accuracy above a tolerance level of five significant digits. There is 
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thus a lower cut-off for the real part of the wavenumber in our results, beyond which 

obtaining numerically accurate results becomes prohibitively expensive computationally. 

The neutral boundary is defined where the imaginary part of the eigenvalue uj and the 

growth rate g are smaller than 1 x 10^^ and 1 x 10^® in magnitude respectively. In 

calculating the streamwise derivatives, nearby stations are taken corresponding to values 

of Reynolds number Re apart by 1% of the value at either of these points. Independence 

of the curve of marginal stability to small deviations in the numerical value of this 

parameter has also been confirmed. 

4.2. Parallel vs. non-parallel analysis 

As mentioned before, the non-parallel approach is expected to deviate significantly from 
the parallel approach in the region where the Reynolds number based on streamwise 
distance is small. 

First, as validation of the OS solver used, we compare the results of the present analysis 
on the h yperbolic tangent profile (th e case of counterflowing streams with A = — 1) with 



those of iBetchov fc Szewczvk figGSl) in figure [3 Secondly, we establish that the result 
showing the flow to be unstable at i?e = is not specific to a counterflow situation, 
rather it is inherent in the parallel flow assumption. Curves of marginal stability from 
a parallel analysis on the similarity proflles for the velocity ratios and 0.5 are also 
plotted in flgureEl These are compared with those from a parallel analysis on a suitably 
transformed hyperbolic tangent profile. 

Note that while the similarity solution has a well-defined length scale (12.1(1 associated 
with it, the hyperbolic tangent profile may be arbitrarily scaled. For the results shown 
in figure 13 the hyperbolic tangent function used is so transformed that its vorticity 
thickness (|4.2(l is identical to that of the similarity solution for the corresponding veloc- 
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Figure 2. Results of OS (parallel) analysis: circles, Betchov & Szewczyk (1963); solid line, 
present analysis on tanh profile; squares, transformed tanh profile with A = 0; short dashes, 
similarity profile with A = 0; crosses, transformed tanh profile with A = 0.5; dash-dot line, 
similarity profile with A = 0.5. 

ity ratio. We note that all the curves approach the origin, irrespective of whether the 
hyperbolic tangent or the similarity profile is used. 

Next, the results of the parallel and non-parallel analyses are compared in figure 13 for 
the four velocity ratios 0.0, 0.25, 0.5 and 0.75. The results correspond to the respective 
similarity profiles. The relevant quantity on the vertical axis for the results of the parallel 
analysis is the physical wavenumber a, whereas for the non-parallel analysis it is the real 
part of the effective wavenumber cioff defined in (|3.()|l . 

Apart from the terms already present in the OS equation, the non-parallel operator 
M includes further terms that essentially provide a correction at low Re. ft is expected 
that the difference between the parallel and non-parallel approaches diminishes in the 
limit Re — + oo. Figure 01 shows that the curves of marginal stability, from both parallel 
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Figure 3. Comparison of the parallel and non-parallel analyses. The Reynolds number is plotted 
on the horizontal axis. The vertical axis corresponds to the real parts of the wavenumber a 
and the effective wavenumber Qoff for the parallel and non-parallel analyses respectively. Short 
dashes, OS; solid line, non-parallel; dash-dot line, Rayleigh. 

and non-parallel analyses, do indeed approach the neutrally stable mode of the solution 
to Rayleigh's equation for each of the velocity ratios considered above. Note that the 
curves for A = have been plotted with Re on a log scale to show that the results are 
indistinguishable as i?e ^ oo; but differences are noticeable even at i?e = 10'^! 

From §131 it is evident that the curve of marginal stability and hence the critical 
Reynolds number can vary depending on the y value at which the growth rate g is 
determined. The monitorin g location makes a huge difference to the stability result in 



the case of boundary layers iJGovindaraian fc Narasimha 1997J) . But for the present flow 



the contribution from the last term in g is smooth in y. The critical Re* ()4.1|) was ob- 
served to vary within a range that was about 10% of the value of the critical Re* at 
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y ^ while traversing from y = —l.O to y ^ 0.5. We do not expect any widely different 

behaviour outside this range as the variation in the velocity profile becomes negligible. 

It is also interesting to note that the minimum of the critical Reynolds numbers for any 

given velocity ratio seems to occur close to the y-location where the streamwise velocity 

profile has the maximum slope (see figure^]) . With y rescaled with respect to the vorticity 

thickness, for the four velocity ratios considered, the location of minimum critical Re* 

deviates most from the location of maximum slope for A = 0.25 (by 0.019 units) and the 

least for A = (by 0.001 units). All the results from the non-parallel analysis presented 

in this paper correspond to a monitoring location fixed at y — 0. 

4.3. Critical Reynolds numbers 

The variation of critical Reynolds number Re with velocity ratio is shown in figure 0] A 
more appropriate Reynolds number in the present flow would be one defined in terms of 
the velocity difference and the vorticity thickness. 

Re* . (1^-^E^ = (i_^i,e (4.1) 

max 

where (1 — X)Uoo is the velocity difference, and 

6Ux)^Ux)^-^-^ (4.2) 

max 

is the vorticity thickness determined by the maximum slope of the velocity profile, '^'max- 

The variation of critical Re* with velocity ratio is shown in figure^ It varies mono- 

tonically with A and approaches a finite limiting value of 55.3 as A goes to 1. Also, the 

variation of critical Re* with A is found to be much less than the variation in critical Re. 
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Figure 4. Variation of the critical Reynolds number with velocity ratio. Circles, 
Re — Uoold/i'', squares, Re* 14.111 . (Dashed and solid curves are fits through the data points.) 

5. Conclusions 

The main result of the present work is the demonstration of the existence of a non- 
zero critical Reynolds number for a laminar plane incompressible mixing layer. Choosing 
the vorticity thickness and velocity differential as length and velocity scales, this critical 
Reynolds number varies in the narrow range from 58.6 at A = to 55.3 at A = 1. The 
present result puts into perspective the prevalent understanding of linear stability of 
the mixing layer. It also underlines the relevance of a non-parallel analysis vis-a-vis a 
parallel one in regard to this flow and other open flows exhibiting a high degree of non- 
parallelism. This is a striking example of a flow where the use of non-parallel theory is 
essential to avoid drawing incorrect physical conclusions. Though parallel flow theory has 
given revealing insights to instability mechanisms for over a century, there are regimes 
of flow where this theory is qualitatively wrong. 
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The finding of a critical Reynolds number which is not too low has an appealing 

consequence. The laminar similarity flow analysed here is physically realizable only at 

some distance downstream of the splitter plate, and the results may therefore be verified 

experimentally or by direct numerical simulations. This is not possible for the earlier 

parallel flow results. We note that since the flow is convectively stable below Re* = 55, 

the question of absolute instability below this Reynolds number does not arise. While 

it is outside the thrust of the present work, it is relevant to mention that fo r inviscid 



flow, from a parallel analysis, it is well established IjHuerre fc Monkewitz 19851) that the 
instability is convective for any mixing layer formed by coflowing streams. 

We also claim that for the purpose of stability analysis, the dependence of results on 
the exact velocity profile is weak, the form suffices to obtain results to a consistently 
good accuracy. For all the non-parallel, parallel and inviscid analyses carried out on the 
similarity profiles, the results differ from corresponding analyses on the (suitably rescaled) 
hyperbolic tangent profiles by no more than 5%. This is in contrast to wall-bounded shear 
flows where the stability results are very sensitive to the mean-flow velocity profile. 

The authors wish to thank the Defence Research and Development Organisation (DRDO), 
India for supporting this work. 
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